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摘 要 : 为 明晰 高 原 流域 土地 利用 变化 时 空 特征 及 其 对 生境 质量 的 影响 ,以 典型 高 原 流域 


拉萨 河流 域 为 例 ， 


基于 GEE 平 台 和 随机 森林 分 类 算法 , 解 译 拉萨 河流 域 1990 一 2020 年 土地 利用 信息 ,分 析 30 a 间 土地 利用 的 时 空 动 


态 变 化 特征 ,并 结合 InVEST 模 型 评估 流域 生境 质量 ,探讨 流域 生境 质量 变化 特征 。 人 研究 结果 表明 :针对 高 寒 高 海 
拨 地 区 的 土地 利用 分 类 问题 ,筛选 植被 生长 季 , 用 相 邻 年 份 的 影像 进行 补充 和 蔡 换 缺失 的 数据 ,并 去 云 后 生成 的 年 
度 合 成 影像 ,提高 了 地 类 分 类 效率 和 精度 ;建立 了 一 种 结合 Landsat 影 像 光谱 特征 .纹理 特征 和 地 形 特 征 的 最 优 RF 
分 类 模型 ,实现 了 对 高 原 流 域 土地 利用 信息 的 提取 ,OA Kappa 系数 和 Fl1-score 分 别 在 88.16% .0.84 和 0.70 以 上 ;20 
世纪 90 年 代 以 来 ,拉萨 河流 域 土地 利用 变化 特征 整体 表现 为 “五 增 三 减 ”, 草 地 ,建设 用 地 林地、 湿地 和 水 域 面积 增 
加 ,耕地 冰川 及 永久 积 雪 和 未 利用 地 面积 减少 ,其 中 ,建设 用 地 的 增长 幅度 最 大 ,达到 288.35% , 而 耕地 缩减 最 多 ， 
达 50.18% , 且 耕 地 减少 主要 是 受到 城镇 范围 扩张 的 挤占 ;未 利用 地 面积 减少 最 为 显著 ,其 转 出 部 分 主要 转化 为 草 
地 ;研究 期 内 高 质量 生境 主要 分 布 在 流域 下 游 .西南 侧 以 及 源头 区 域 ,以 草地 水体. 冰川 及 永久 积 雪 和 湿地 等 土地 


= 


re 


j 类 型 为 主 ; 而 拉萨 市 市 辖区 , 达 孜 区 PK] EL Hk Ee TREHI , DA Betta EDR AS A SE 


域 生境 质量 较 差 ;流域 整体 生境 质量 呈 "* 先 减 后 增 再 减 然后 趋 缓 " 的 变化 趋势 ,生境 质量 指数 从 0.53 提升 到 0.57; 人 


口 增加 和 快速 城镇 化 加 快 了 城镇 范围 扩张 和 耕地 的 减少 ,但 气候 变化 和 生态 治理 工程 对 拉萨 河流 域 的 生境 质量 改 


善 均 起 到 了 积极 作用 ;本 研究 可 为 制定 高 原 地 区 流域 生态 修复 策略 .促进 土地 利用 可 持续 以 及 提升 生境 质量 提供 


科学 依据 与 参考 。 
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四 由 


生境 质量 是 衡量 区 域 生态 环境 为 个 体 或 种 群 
持续 发 展 提供 生存 条 件 能 力 的 衡量 标准 ,是 区 域 生 
物 多 样 性 水 平 的 具体 反映 。 但 是 ,全 球 气候 变化 
和 人 类 活动 之 来 的 土地 利用 变化 严重 影响 区 域 生 
境 的 时 空 分 布 格局 ,从 而 对 生境 质量 产生 影响 ”。 
因此 ,通过 解 译 不 同 空 间 尺 度 的 土地 利用 信息 并 评 
估 其 生境 质量 ,进而 揭示 区 域 土地 利用 和 生境 质量 
变化 规律 及 空间 分 异 特征 ,对 生态 文明 建设 下 的 空 
间 治 理 具有 重要 意义 。 

青藏 高 原作 为 对 全 球 气候 和 人 类 活动 变化 响 
应 最 敏感 的 区 域 ”, 具 有 生态 环境 脆弱 .经 济 社会 发 
展 较 落后 等 特征 ”, 在 高 原生 态 修复 与 土地 整治 研 
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究 方 面具 有 典型 性 。 已 有 人 研究 表明 ,青藏 高 原 土 地 
利用 空间 布局 差异 显著 ,人 口 随 流域 分 布 且 相 对 集 
中 中 ,流域 作为 自然 要 素 与 人 文 要 素 相 互 关联 相互 
作用 的 地 理 单元 和 城市 发 展 过 程 中 人 类 活动 的 载 
体 ” ,具有 良好 的 生态 完整 性 ”, 但 社会 经 济 快速 发 
展 和 快速 城镇 化 为 其 带 来 了 日 益 严 重 的 环境 问 
题 。 土 地 利用 格局 不 断 变化 ,生态 环境 压力 不 断 升 
级 ,导致 流域 生态 系统 产生 景观 空间 破碎 ”土地 退 
化 水 资源 短缺 路 和 水 环境 污染 "等 问题 。 探 究 
流域 土地 利用 变化 ,提升 流域 生境 质量 ,有 利于 优 
化 流域 土地 利用 开发 保护 格局 ,确保 生态 系统 健康 
可 持续 。 而 有 关 青 藏 高 原 土 地 利用 和 生境 质量 变 
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化 的 研究 , 仅 有 部 分 学 者 采用 大 尺度 AIREY) 
土地 利用 数据 进行 定量 描述 ”或 者 仅 针 对 某 单 一 
土地 利用 类 型 的 时 空 变 化 进行 研究 ” ,忽视 了 城市 
流域 具有 土地 利用 多 样 性 和 生态 系统 完整 性 在 该 
研究 领域 的 典型 性 。 

拉 椭 河流 域 地 处 青藏 高 原 中 南部 ,位 于 西藏 自 
治 区 社会 .经济 文化 中 心地 带 , 是 典型 的 高 原 城市 
流域 。 本 人 研究 基于 GEE(Google Earth Engine) 平 台 ， 
选取 拉萨 河 流域 1990 一 2020 年 7 期 Landsat 系列 影 
像 数 据 ,构建 光谱 纹理 和 地 形 3 种 分 类 特征 ,利用 
随机 森林 (Random Forest, RF ) 算 法 进行 土地 利用 分 
类 ,通过 分 析 流 域 30a 长 时 间 序 列 下 的 土地 利用 变 
化 特征 ,结合 InVEST 模 型 评估 研究 期 内 生境 质量 变 
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化 , 旨 在 探索 RF 方 法 和 Landsat 数 据 在 高 原 城 市 流 
域 土 地 利用 分 类 中 应 用 的 可 行 性 ,为 制定 青藏 高 原 
流域 生态 恢复 策略 ,提升 流域 生态 系统 服务 功能 ， 
促进 土地 利用 可 持续 发 展 提供 科学 依据 。 


1 数据 与 方法 


1.1 研究 区 概况 

LB EIA fit dk (90° 05'~93°20'E, 29°20'~31°15'N) 
位 于 西藏 自治 区 中 南部 (图 1) ,是 雅鲁藏布江 的 一 
级 支流 ,流域 面积 32548 km’, 属 高 原 温带 半 干 旱季 
风气 候 , 多 年 平均 降雨 量 约 500 mm ,多 年 平均 气温 
约 为 7 °C ,平均 相对 湿度 约 45% , 流域 北部 地 形 平 
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图 1 拉萨 河 流域 概况 及 高 分 辩 率 影像 验证 样本 点 分 布 


Fig. 1 Overview of Lhasa River Basin and distribution of high-resolution image verification sample points 
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坦 , 往 南 地 形 起 伏 度 逐 渐 增 大 ,包含 山地 ` 盆 地 和 河 
谷 平原 等 地 形 ,流域 内 分 布 有 季节 性 冻 土 和 冰川 ， 
最 大 积 雪 厚 度 达 11 cm, 植 被 类 型 主要 是 高 山 草 馈 
和 高 山 草原 等 ,垂直 分 带 较为 明显 ,是 西藏 自治 区 
人 口 密集 区 和 经 济 发 展 的 核心 区 以 及 重要 的 农 产 
品 生产 区 。 
1.2 数据 来 源 与 处 理 
1.2.1 数据 来 源 与 获取 GEE 是 一 个 基于 网 络 的 地 
理 空间 分 析 和 数据 处 理 平台 ,储存 了 PB 级 的 遥感 影 
像 和 地 球 科学 数据 集 ,为 科研 人 员 提 供 了 一 个 具有 
行星 级 运算 处 理 能 力 和 内 在 并 行 处 理 方式 的 数据 
采集 、 人 处 理 和 分 析 的 平台 ""“"。 本 研究 使 用 Landsat 
TM、OLI 地 表 反 射 率 数据 集 和 SRTM 数字 高 程 数 据 
来 解 译 拉萨 河流 域 1990 年 1995 年 .2000 年 .2005 
年 .2010 年 .2015 年 和 2020 年 7a 的 土地 利用 数据 ， 
所 有 数据 的 空间 分 辨 率 均 为 30 mo 

Landsat TM 、OLI 地 表 反 射 率 数 据 和 SRTM 数字 
高 程 数 据 的 所 有 处 理 均 在 GEE 平 台 (https://earthen- 
gine.google.com) 上 进行 。 主 要 包括 以 下 几 个 步骤 : 
(1) 为 了 纠正 高 寒 高 海拔 地 区 光学 影像 多 云 的 问 
题 , 选 取 了 每 年 植被 生长 季 (5 一 9 月 ) 的 地 表 反 射 率 
数据 ,用 每 个 研究 年 份 前 后 的 影像 奉 换 或 补充 云雾 
覆盖 的 影像 ,分 别 得 到 了 1990 年 .1995 年 .2000 年 、 
2005 年 的 77 景 .143 景 .133 景 .216 景 Landsat 5 TM 
影像 和 2010 年 、2015 年 271 景 、340 景 Landsat 7 
ETM+ 影 像 和 2020 年 413 景 Landsat 8 OLI 影 像 。(2) 
对 云 量 小 于 20% 的 影像 进行 了 去 云 处 理 。(3) 利用 
reducer 也 数 对 影像 进行 中 值 合成 ,得 到 每 个 研究 年 
的 单 幅 影 像 .(4) 构建 分 类 特征 集 来 提高 分 类 精度 ， 
计算 每 幅 影 像 的 归 一 化 植被 指数 (NDVI)'" 1、 归 一 化 
建筑 指数 (NDBI)'\ 归 一 化 水 指数 (NDWI)"" 植被 


水 分 含量 指数 (LSWI)"”"。 利 用 灰 度 共生 矩阵 
(Gray-level Co-occurrence Matrix, GLCM ) 提取 影像 
纹理 特征 ,用 glemTexture KAOTE T fA ZIE Xt 
比 度 .相关 性 方差 EE IE 6AA I CE 
MEU ;并 提取 了 STRMGCL1_003 数字 高 程 数据 的 海 
R .坡度 、 坡 向 、 山 体 阴影 等 4 个 地 形 特征 。 上 述 步 
又 完成 后 ,得 到 了 每 个 研究 年 用 于 土地 利用 解 译 的 
IH RY RA 
12.2 分 类 与 精度 评价 ”参照 (土地 利用 现状 分 类 》 
(GB/T21010-2017) ,结合 研究 区 土地 利用 现状 和 特 
点 ,将 土地 利用 类 型 分 为 水 体 、 林 地 草地、 耕地 、 建 
设 用 地 ,未 利用 地 .湿地 和 冰川 及 永久 积 雪 8 类 。 本 
人 研究 采用 融合 多 特征 要 素 的 土地 利用 分 类 框架 (图 
2) ,根据 不 同 地 物 在 不 同 波段 组 合 下 存在 的 特征 差 
异 (包括 形状 .颜色 纹理 ) ,采用 在 GEE 上 直接 标注 
的 方法 选取 用 于 土地 利用 分 类 的 样本 数据 , 共 选 取 
样本 总 数 17472 个 (其 中 80% 作 为 训练 样本 ,20% 用 
于 精度 验证 )。 为 了 进一步 验证 地 物 分 类 的 准确 
性 ,本 研究 使 用 随机 采样 方法 在 Google Earth 历史 高 
分 辨 紊 影像 中 收集 了 所 有 人 研究 年 份 共 611 个 随机 验 
证 点 对 结果 进行 评估 (图 1)。 

研究 表明 ,相同 条 件 下 ,RF 随机 森林 算法 对 数 
据 品 声 和 过 拟 合 有 较 好 的 处 理 能 力 ,并 且 相 较 于 
其 他 机 器 学 习 算法 的 分 类 精度 较 高 ” ,因此 ,本 研 
究 选择 了 RF 作为 分 类 带 。 决 策 树 数量 对 分 类 精度 
有 具 有 重要 的 影响 ” ,综合 考虑 分 类 精度 和 硬件 性 
能 ,将 决策 树 数 量 设置 为 ntree=50, 基于 获取 的 样 
本 ,在 GEE 平 台中 使 用 ee. Classifier. smileRandomFor- 
est 函数 ,获得 每 个 研究 年 的 土地 利用 分 类 图 。 
1.3 研究 方法 
1.3.1 土地 利用 动态 度 ”土地 利用 动态 度 分 为 单一 
土地 利用 动态 度 (L) 和 综合 土地 利用 动态 度 (L.)。 
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>) 土地 利用 动态 度 、 土 地 利用 程度 、 土 地 利用 转移 矩阵 30 a 间 土地 利用 变化 
土地 利用 图 谱 ( InVEST 生 境 质量 模块 => 30 a 间 生境 质量 变化 
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图 2 技术 路 线 


Fig. 2 Technical route 
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单一 土地 利用 动态 度 用 来 表征 研究 区 内 某 土地 利 
用 类 型 在 一 定时 间 范围 内 的 数量 变化 情况 ;综合 二 
地 利用 动态 度 则 表征 在 一 定时 间 范围 内 研究 区 所 
有 土地 利用 类 型 的 整体 数量 变化 情况 ,其 计算 公式 
分 别 为 = 


p= tax 1x 100% (1) 
SaL, 
[= 人 一直 x 100% (2) 


式 中 : 工 为 单一 土地 利用 动态 度 ;UV, 和 0 分 别 为 研究 
区 初期 和 末期 某 土地 利用 类 型 的 面积 ; 埃 是 综合 十 
地 利用 动态 度 ;L 是 研究 初期 第 i 类 土地 利用 类 型 
的 面积 ; AL 是 研究 期 内 第 i 类 土地 利用 类 型 转变 
为 非 i 类 土地 利用 类 型 面积 的 绝对 值 ;7 为 研究 期 
(年 )。 
1.3.2 土地 利用 程度 土地 利用 程度 综合 反应 了 在 
自然 因素 和 社会 因素 的 共同 影响 下 ,对 土地 资源 利 
用 的 广度 和 深度 。 土 地 利用 程度 综合 指数 依据 不 
同 地 类 进行 分 级 赋值 (本 研究 将 未 利用 地 、 冰 川 及 
KARE 湿地 赋值 为 1; 林 地 草地 水体 为 2; 耕 地 
为 3; 建 设 用 地 为 4) ,用 来 表征 研究 区 土地 利用 的 集 
约 程度 ,其 计算 公式 为 

L=100x YAxC, Lel100,400] (3) 


式 中 :为 研究 区 土地 利用 程度 综合 指数 ;4 为 研究 
区 内 第 i 级 土地 利用 程度 分 级 指数 ;C; 为 研究 区 内 第 
i 级 土地 利用 程度 分 级 面积 百分比 ;n 为 土地 利用 程 
度 分 级 数 。 

1.3.3 土地 利用 转移 矩阵 ”土地 利用 转移 矩阵 可 用 
来 刻画 和 描述 区 域 土地 利用 变化 的 结构 特征 和 土 
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地 利用 类 型 变化 的 方向 ,其 数学 形式 为 : 


Si, Sy Sis S, 
Sx Sy Sa S, 

S;= Sa Sy) Sa; na San (4) 
Sat Six Sia me Sa 


式 中 :5 代表 面积 ;n 代 表土 地 利用 类 型 数 ;iy 分 别 代 
表 研 究 初 期 与 末期 的 土地 利用 类 型 。 

1.3.3 InVEST 生境 质量 模块 ”根据 不 同 生 境 对 胁 
迫 因 子 的 敏感 度 不 同 , 建 立 不 同 地 类 与 胁迫 因子 间 
的 联系 ,通过 计算 生境 质量 指数 来 评估 研究 区 在 每 
个 研究 年 的 生境 质量 空间 分 布 及 变化 情况 , 取 值 范 
于 介 于 0~1 ,数值 越 大 ,代表 生境 质量 越 高 所。 模型 
中 涉及 的 主要 参数 包括 土地 利用 数据 .胁迫 因子 数 
据 、 各 胁迫 因子 的 最 大 影响 距离 和 距离 衰减 函数 、 
各 胁迫 因子 的 权重 、 不 同 土地 覆盖 类 型 的 生境 适宜 
度 及 其 对 胁迫 因子 的 敏感 度 。 参 考 相关 文献 [25- 
27] ,并 根据 拉萨 河流 域 实际 情况 及 专家 建议 进行 
相应 参数 赋值 ( 表 1, 表 2)。 


2 结果 与 分 析 


2.1 分 类 精度 评价 
分 类 结果 表明 ,1990 年 1995 年 .2000 年 .2005 
Æ 2010 年 .2015 年 和 2020 年 总 体 精 度 分 别 为 


R1 胁迫 因子 属性 表 
Tab.1 Attribute table of stress factors 


胁迫 因子 。 最 大 影响 距离 /km? 权重 衰减 函数 
耕地 6 0.7 线性 
建设 用 地 10 1 指数 
未 利用 地 5 0.3 Ha 


R2 不 同 土地 利用 类 型 的 生境 适宜 度 及 其 对 胁迫 因子 的 敏感 度 
Tab. 2 Habitat suitability of different land use types and their sensitivity to stress factors 


对 胁迫 因子 的 敏感 度 


cee 耕地 建设 用 地 未 利用 地 
冰川 及 永久 积 雪 1 0.6 0.8 0.3 
草地 0.79 0.3 0.3 0.6 
耕地 0.25 0.3 0.3 0.4 
建设 用 地 0 0 0 0 
林地 0.68 0.8 0.9 0.5 
湿地 1 0.7 0.9 0.4 
水 体 1 0.8 0.8 0.3 
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90.5% , 92.18% .88.16% , 91.22% 92.72% 、96.19% 、 
93.78% , Kappa 系数 分 别 为 0.88、0.91、0.84、0.89、 
0.91 .0.95 .0.92 , 混 消 矩 阵 显 示 每 种 土地 利用 类 型 的 
详细 分 类 精度 ( 表 3) ,总 体 分 类 精度 满足 土地 利用 
分 类 的 要 求 。 此 外 ,为 了 进一步 验证 分 类 精度 ,将 
其 与 高 分 辩 率 影像 上 选取 的 样本 进行 了 比较 , 结 
表明 :对 冰川 及 永久 积 雪 草地、 林地 和 湿地 的 分 类 
准确 度 较 好 (Fl-score 分 别 为 0.96, 0.91, 0.92 和 
0.91) ,耕地 和 水 体 次 之 (F1-score 分 别 为 0.82 和 
0.88) ,未 利用 地 和 建设 用 地 稍 差 (Fl1-score 分 别 为 
0.72 和 0.70)。 

2.2 拉萨 河流 域 土地 利用 变化 分 析 

2.2.1 土地 利用 时 空 动态 拉萨 河流 域 土地 利用 类 
型 以 草地 和 未 利用 地 为 主 , 研 究 期 各 地 类 变化 特征 


具体 表现 为 "五 增 三 减 ”, 草 地、 建设 用 地 、 林 地 ` 湿 
地 和 水 域 面积 增加 ,耕地 、 冰 川 及 永久 积 雪 和 未 利 
用 地 面积 减少 (图 3, 图 4)。 其 中 ,草地 作为 拉萨 河 
流域 面积 最 大 的 地 类 ,1990 一 2020 年 面积 占 流 域 总 
面积 的 比例 均 高 于 49% ,广泛 分 布 于 流域 内 ;从 面 
积 变化 幅度 而 言 ,建筑 用 地 面积 增长 幅度 最 大 ,从 
1990 年 的 22.57 km ,增长 至 2020 年 的 87.65 km?, 增 
长 幅度 高 达 288.35% ,集中 分 布 在 拉萨 市 辖区 ;耕地 
面积 缩减 最 快 , 达 -50.18% , 主要 分 布 在 流域 中 下 游 
的 林 周 县 南部 地 区 .黑人 竹 工 卡 县 中 西部 AAR \ 城 
关 区 和 曲 水 县 。 

2.2.2 土地 利用 动态 度 ”人 研究 期 流域 土地 利用 综合 
动态 度 呈 “ 先 增 后 减 ” 的 变化 趋势 ,在 2005 一 2010 年 
达到 最 大 值 3.56% ,表明 该 时 期 人 类 活动 对 土地 利 


表 3 分 类 精度 


Tab. 3 Classification accuracy 


ah, 19904F 1995 4F 2000 年 2005 年 2010 4F 2015 4F 20204 
土地 利用 类 型 
PA UA PA UA PA UA PA UA PA UA PA UA PA UA 
水 体 91.89 97.14 90.63 93.55 89.19 94.29 95.24 95.24 92.00 97.87 97.67 95.45 97.14 97.14 
林地 93.02 95.24 97.62 95.35 92.96 88.00 95.83 94.52 93.64 94.19 96.98 98.09 96.80 93.81 
草地 90.36 85.23 89.71 91.04 90.87 86.70 91.71 84.69 93.68 90.20 98.15 95.13 94.64 92.71 
耕地 91.11 89.13 88.89 90.24 78.23 94.87 80.77 93.33 88.37 90.48 89.87 98.61 84.38 94.74 
建设 用 地 75.00 82.46 90.00 69.23 82.35 93.33 92.00 95.83 90.00 96.43 90.91 78.96 87.88 96.67 
未 利用 地 96.34 89.77 9848 92.86 91.67 86.43 93.96 93.44 95.53 94.38 97.27 95.03 95.31 95.81 
冰川 及 永久 积 雪 66.67 97.68 71.43 89.54 73.84 87.50 86.67 92.75 66.67 76.92 8545 9742 86.15 93.45 
湿地 79.31 92.00 83.33 89.28 70.83 89.47 78.13 87.68 81.08 93.75 94.74 92.31 90.00 88.52 
OA 90.50 92.78 88.16 91.22 92.72 96.19 93.78 
Kappa 0.884 0.909 0.844 0.886 0.905 0.949 0.918 
TE PA 为 生产 者 精度 ;UA 为 用 户 精 度 ;OA 为 总 体 精度 。 
(a) 1990 年 (b) 1995 年 (c) 2000 年 土地 利用 类 型 
冰川 及 永久 积 雪 
N Mm He, 
A MERA 
E KH 
EE ki% 
E e 
50 km 50 km 0 50km ”上 四 耕地 
国 四 草地 


(d) 2005 年 (e) 2010 年 


50 km 0 


| 


50 km 


(f) 2015 年 (g) 2020 年 


50 km 50 km 


图 3 1990 一 2020 年 拉萨 河流 域 土地 利用 图 谱 
Fig.3 Land use map of Lhasa River Basin from 1990 to 2020 
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图 4 1990 一 2020 年 拉萨 河流 域 土地 利用 结构 
Fig.4 Land use structure in Lhasa River Basin from 
1990 to 2020 


用 的 影响 最 为 强烈 ( 表 4)。 从 单一 土地 利用 动态 度 
来 看 ,30 a H] ,建设 用 地 面积 持续 增加 ,1990 一 1995 
年 增幅 最 大 ,达到 22.11% ,随后 快速 减 小 ,2015 一 
2020 年 仅 为 0.1% ,表明 城镇 化 由 快速 转 为 平稳 推 
进 的 状态 。 人 研究 期 内 ,草地 面积 变化 呈 “ 减 一 增 
增 一 减 一 增 一 减 ”, 与 林地 和 未 利用 地 面积 变化 趋 
势 相 对 应 ,表明 它们 之 间 的 转化 可 能 具有 一 致 性 。 
冰川 及 永久 积 雪 在 2000 一 2005 年 面积 大 幅 增加 , 达 
到 45.58% ,相关 研究 所 表明 , 2000—2004 年 拉萨 河 
流域 气温 出 现 较 大 幅度 降低 ,从 而 导致 冰川 及 永久 
积 雪 面 积 增加 。 湿 地 和 水 体面 积 总 体 呈 增加 趋势 ， 
研究 期 内 分 别 增加 40.4% 和 12.83% ,全 球 气候 变 暖 
加 速 冰雪 消融 是 该 地 类 面积 增加 的 主要 原因 。 

2.2.3 土地 利用 程度 流域 土地 利用 程度 综合 指数 
是 “ 先 增 后 减 ” 变 化 ,其 值 介 于 162.20~171.15 ,2005 
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年 达到 最 大 值 ,随后 有 所 减 小 (图 $)。 土 地 利用 程 
度 随 海拔 高 度 的 增加 呈 先 下 降 后 趋 缓 特 征 ,城关 
区 JAK . 堆 龙 德 庆 区 ` 林 周 县 . 墨 竹 工 卡 县 . 曲 水 
县 等 土地 利用 程度 相对 较 高 ,与 这 些 区 域 平 均 海 拔 
相对 较 低 ( 介 于 3500~4500 m) 有 关 , 其 人 类 经 济 建 
设 生 产 活动 频繁 、 土 地 利用 综合 开发 程度 相对 较 
高 ; 当 雄 县 RE ` 色 尼 区 等 高 海拔 区 域 土地 利用 
程度 较 低 , 分 别 为 157.26、159.21、152.31; 平 均 土 地 
利用 程度 综合 指数 均 低 于 160, 而 比如 县 、 桑 日 县 的 
部 分 区 域 基本 为 草地 ,导致 土地 利用 程度 综合 指数 
偏 高 。 

2.2.4 土地 利用 转移 特征 ”人 研究 期 内 流域 各 地 类 之 
间 转 换 频 繁 (图 6, 表 5)。 其 中 未 利用 地 是 草地 转 入 
面积 最 多 的 地 类 ,说 明 流 域 近年 来 生态 修复 与 绿化 
工程 发 挥 了 重要 作用 ;建设 用 地 转 人 面积 远大 于 转 
出 面积 ,主要 是 为 满足 人 口 向 城镇 集聚 .生产 生活 
用 地 需求 增加 ,大 面积 的 草地 、 耕 地 和 未 利用 地 被 
挤占 与 开发 ,三 者 分 别 向 建设 用 地 转 人 了 28.96 
km? 19.46 km? 和 19.38 km’; 冰川 及 永久 积 雪 转 出 面 
积 远大 于 转 人 面积 ,再 次 表明 全 球 气 候 变 暖 加剧 了 
高 原 区 域 的 冰雪 消融 ;林地 转 人 面积 约 为 转 出 面积 
的 2 倍 ,说 明 流域 多 年 实施 的 退耕 还 林 还 草 工 程 取 
得 了 显著 成 效 ;湿地 和 水 域 转 和 人 面积 远大 于 转 出 面 
只 ,说 明 流域 内 水 环境 治理 得 到 了 较 好 的 效果 ; 综 
合 表 明 流 域 生 态 环境 治理 效果 较 好 。 

2.2 拉萨 河流 域 生境 质量 变化 

2.2.1 生境 质量 空间 分 布 根据 模型 计算 结果 (图 
7),1990 年 .1995 年 .2000 年 .2005 年 .2010 年 .2015 
年 和 2020 年 拉萨 河流 域 的 生境 质量 指数 平均 值 分 
别 为 0.53、0.51、0.56、0.61、0.56、0.58 Fil 0.57 ,整体 呈 


表 4 不 同时 期 拉萨 河流 域 土地 利用 动态 度 指数 


Tab.4 Land use dynamic index of Lhasa River Basin in different periods 1% 
土地 利用 类 型 1990 一 1995 年 1995 一 2000 年 。 2000 一 2005 年 。 2005 一 2010 年 。 2010 一 2015 年 。 2015—20204F 

冰川 及 永久 积 雪 -7.03 -5.83 45.58 -11.38 -1.02 -3.53 
草地 -0.74 4.43 0.95 -2.26 0.72 -0.82 
耕地 -6.13 -6.07 -1.39 3.84 1.47 -2.68 
建设 用 地 22.11 4.71 4.06 3.39 1.12 0.1 

林地 1.94 -7.75 -4.14 17.85 -0.54 4.4 

湿地 4.23 -8.28 -0.57 12.55 2.92 1.83 
水 体 0.03 3.59 -2.34 -4.86 1.55 6.53 
未 利用 地 1.53 -3.29 -3.84 4.89 -1.53 0.52 


综合 土地 利用 动态 度 2.95 3.35 


3.35 3.56 1.99 2.00 
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图 5 全 流域 及 各 区 (县 ) 土 地 利用 程度 综合 指数 


Fig. 5 Comprehensive index of land use degree of the whole basin and districts (counties) 
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图 6 1990 一 2020 年 拉萨 河流 域 土地 利用 变化 转移 方向 
Fig.6 Transfer direction of land use change in Lhasa River Basin from 1990 to 2020 


“ 先 减 后 增 再 减 然 后 趋 缓 " 变 化。 流域 下 游 .西南 侧 
以 及 源头 区 域 生境 质量 较 好 ,以 草地 、 湿 地 和 冰川 
及 永久 积 雪 等 土地 利用 类 型 为 主 ; 拉 蔷 市 市 辖区 、 
BAK 、 林 周 县 . 曲 水 县 和 墨 竹 工 卡 县 县 城 周边 ,以 
及 流域 中 上 游 等 区 域 生 境 质量 较 差 ;1990 一 2000 年 
生境 质量 空间 格局 基本 未 发 生变 化 ,而 2000 一 2005 
年 流域 中 部 生境 质量 大 面积 提升 ,原因 是 中 部 冰川 


及 永久 积 雪 面积 大 幅 增 加 ;2010 一 2020 年 ,流域 生 
境 质量 空间 分 布 格局 逐渐 恢复 且 生 境 质 量 明显 
提升 。 

2.2.2 生境 质量 转移 特征 根据 生境 质量 评估 结 

R ,将 拉萨 河流 域 生境 质 量 划 分 为 差 (0~0.2) 一般 
(0.2~0.5) 、 良 好 (0.5~0.8) 和 优 (0.8~1) 共 4 个 级 别 。 
1990 一 2020 年 流域 11.72% 的 区 域 生 境 质 量 等 级 升 
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表 5 1990 一 2020 年 拉萨 河流 域 土地 利用 转移 矩阵 
Tab.S Land use transfer matrix of Lhasa River Basin from 1990 to 2020 


2020 4E 
Ga 冰川 及 永久 积 雪 草地 耕地 ”建设 用 地 ”林地 ”湿地 ”水 休 ”未 利用 地 总 计 
冰川 及 永久 积 雪 - 85.06 0 0 0.00 1.60 0.63 604.68 691.98 
草地 45.99 - 50.24 28.96 762.01 433.41 49.78 420.29 1790.68 
耕地 0 306.64 - 19.46 4.62 7.27 13.88 20.11 371.98 
建设 用 地 0 2.92 0.93 一 0.03 0.14 1.44 1.09 6.55 
林地 0 336.67 0.16 0.46 一 24.74 0.81 73.86 436.70 
湿地 0 176.91 0.44 0.21 14.84 - 9.15 1.85 203.41 
水 体 2.46 32.54 3.12 3.17 1.33 14.72 一 38.03 95.36 
未 利用 地 47.23 2301.56 33.38 19.38 173.85 5.80 54.36 一 2635.55 
总 计 95.68 3242.30 88.28 71.63 956.69 487.67 130.05 1159.91 - 
(a) 19904 (b) 19954 (c) 20004 生境 质量 指数 
N 002 
人 HE 0.2~0.5 
0.5~0.8 
E 0.8~1 
0 50km 0 50km 0 50km 
(d) 2005 年 (e) 2010 年 (f) 20154 (g) 2020 年 


0 50km 


图 7 1990 一 2020 年 拉萨 河流 域 生 境 质量 空间 分 布 及 其 变化 
Fig.7 Spatial distribution and change of habitat quality in Lhasa River Basin from 1990 to 2020 


高 ,向 “良好 ”和 “ 优 ” 转 换 ,3.93% 的 区 域 生境 质量 下 
降 , 以 “良好 ”和 “ 优 ” 向 “ 差 " 转 换 为 主 ( 表 6);1990 一 
2000 年 8.02% 的 生境 状况 得 到 改善 ( 差 和 一 般 向 良 
好 和 优 转 换 ) ,4.04% 的 区 域 生境 质量 水 平 降低 ; 
2000 一 2010 年 转换 为 “良好 ”和 “ 优 ” 的 面积 占 比 分 
别 为 2.72% 和 0.96% , 低 等 生境 质量 恢复 缓慢 ,与 第 
一 个 10a 的 生境 状况 基本 保持 不 变 , 只 有 小 部 分 生 
境 质 量 上 升 ;2010 一 2020 年 4.4% 的 区 域 生 境 质量 水 
平 提高 ,而 2.86% 的 区 域 下 降 ;整体 而 言 ,1990 一 
2020 年 流域 的 生境 状况 有 所 提升 ,质量 下 降 区 主要 
位 于 城市 及 其 周边 , 拉 艾 市 市 边区 下 降 最 为 显著 是 
呈现 明显 的 聚集 状态 ,可 能 是 由 于 人 口 增长 ,城市 
建成 区 扩张 ,建设 用 地 面积 增加 ,进而 对 周围 生境 
质量 造成 威胁 ,致使 该 区 域 生境 质量 下 降 。30 a 间 
流域 内 有 19.15% 的 区 域 土地 利用 类 型 发 生 了 变化 ， 
主要 发 生 在 未 利用 地 向 草地 的 转移 , 相 较 于 1990 年 
2020 年 草地 面积 增加 了 1475.64 km?。 建 设 用 地 、 耕 


地 等 生境 胁迫 因子 未 发 生 较 大 变化 ,建设 用 地 H 
地 变化 与 生境 质量 变化 的 最 大 相关 性 仅 为 0.001、 
0.003 ,因此 ,研究 期 内 生境 质量 发 生变 化 的 主导 因 
子 是 未 利用 地 面积 的 减少 ,流域 整体 生境 质量 也 从 
0.53 提升 到 0.57。 


3 讨论 


如 前 文 所 述 ,2 种 分 类 精度 验证 方法 均 显示 本 
研究 生产 的 拉萨 河流 域 土地 利用 数据 准确 可 靠 , 构 
建 的 14 种 分 类 特征 、 采 用 Landsat 系列 数据 和 RF 算 
法 在 高 寒 高 海拔 地 区 的 土地 利用 分 类 中 的 应 用 是 
可 行 的 ,与 伐 蒙 京 等 ”对 高 寒 湿地 地 区 土地 覆盖 分 
类 方法 的 研究 结果 是 一 致 的 。 但 是 仍然 存在 一 定 
的 局 限 性 ,一 是 生产 的 土地 利用 数据 的 质量 受到 
Landsat 影像 质量 和 数量 的 影响 ,由 于 高 海拔 地 区 光 
照 强 度 较 高 ,沙化 地 和 不 透水 面 旦 现 相似 的 反射 特 
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表 6 1990 一 2020 年 生境 质量 等 级 转移 矩阵 
Tab. 6 Habitat quality grade transfer matrix from 1990 to 2020 

IN BE 生境 质量 等 级 差 一 般 良好 优 
1990 一 2000 年 差 - 10.57 2166.61 125.04 
一 般 11.75 一 295.88 21.58 
良好 587.75 25.42 一 234.86 

优 699.25 1.06 491.33 一 
2000 一 2010 年 差 一 18.10 835.27 310.36 
一 般 21.47 一 51.32 3.19 
良好 968.82 76.38 一 475.97 

优 248.65 11.30 115.19 - 
2010—2020 4 差 - 27.75 1132.69 197.46 
一 般 18.71 一 98.06 4.95 
良好 564.91 65.17 一 394.84 

优 292.00 7.56 179.55 - 
1990—20204F 差 - 34.40 2497.44 110.28 
一 般 39.76 一 311.14 21.26 
良好 543.06 50.57 - 559.03 

优 648.76 3.59 315.85 - 


征 ,导致 本 研究 未 利用 地 和 建设 用 地 分 类 时 易 误 分 
或 错 分 ;其 次 ,本 研究 中 地 物 分 类 使 用 的 训练 集 和 
验证 集 均 为 人 工 目 视 采集 ,会 给 数据 生产 和 精度 评 
佑 带 来 一 定 程度 的 不 确定 性 。 

土地 利用 变化 与 人 类 活动 密切 相关 ,生活 和 生 
产 增加 对 不 同 土地 利用 类 型 的 需求 ,人 口 增长 和 城 
市 化 增加 建设 用 地 转 入 面积 ,开发 建设 活动 驱动 有 
限 的 土地 资源 频繁 转换 。1990 一 2020 年 拉萨 市 常 
住人 口 增加 51.06x10* 人 ;2006 一 2020 年 GDP 增 长 
575.77x10: 元 ,旅游 业 的 快速 发 展 、 二 三 产业 比例 不 
斯 攀升 ,导致 城镇 范围 扩张 ,挤占 周围 耕地 ;30 a 间 
草地 和 未 利用 地 之 间 的 相互 转换 面积 达 2721.85 
km ,是 变化 量 最 大 的 地 类 转换 ,与 陈 伊 多 等 对 西 
藏 自治 区 土地 利用 变化 的 研究 结果 一 致 。1995 一 
2018 年 青藏 高 原 草 地 和 建设 用 地 持续 增加 ,耕地 、 
未 利用 地 及 冰川 面积 减少 “ ,于 本 研究 得 出 的 
1990 一 2020 年 拉萨 河流 域 土地 利用 变化 表现 为 “五 
增 三 减 ” ,草地 建设 用 地 林地 .湿地 和 水 域 面积 增 
加 ,耕地 、 冰 川 及 永久 积 雪 和 未 利用 地 面积 减少 结 
果 相 近 。 此 外 ,本 研究 表明 ,1990 一 2020 年 拉萨 河 
流域 整体 生境 质量 由 0.53 HEFT 28 0.57 , 2005 年 达到 
最 大 值 ,2010 年 .2015 年 和 2020 年 整体 生境 质量 分 
别 为 0.56.0.358 和 0.57, 较 2005 年 略 有 下 降 , 与 税 燕 
PESO 对 生境 质量 的 研究 结果 有 偏差 ,这 可 能 是 由 
于 空间 分 辩 率 .数据 源 和 分 类 方法 的 不 同 导致 不 同 


的 土地 利用 产品 之 间 存 在 较 大 差异 造成 的 。 

目前 ,对 生境 质量 的 评估 一 般 从 两 个 方面 进行 
考虑 ,一 是 土地 利用 类 型 自身 的 生境 质量 差异 ,二 
是 人 类 活动 等 对 生境 质量 的 威胁 。 且 多 数 人 研究 基 
于 InVEST 模 型 进行 生境 质量 评估 ,就 模型 可 徘 性 而 
言 ,存在 一 些 不 足 :(1) 模型 将 胁迫 因子 的 威胁 县 
加 ,造成 多 种 胁迫 因子 的 威胁 远 高 于 个 体 威胁 之 
和 ,计算 的 生境 质量 存在 误差 ;(2) 人 研究 区 的 生境 质 
量 一 般 是 受 内 外 两 部 分 的 生境 胁迫 ,而 在 InVEST 模 
型 中 往往 忽略 了 人 研究 区 范围 之 外 的 生境 威胁 ,计算 
结果 可 能 存在 威胁 强度 降低 的 现象 。 

本 人 研究 以 拉萨 河流 域 为 研究 对 象 , 基 于 CEE 平 
台 ,利用 多 源 遥 感 数据 和 机 咒 学 习 分 类 算法 对 流域 
土地 利用 类 型 进行 解 译 ,并 用 解 译 结果 评估 了 流域 
生境 质量 ,最 终 取得 了 一 定 研 究 成 果 ,结果 可 为 流 
域 土 地 可 持续 利用 制定 生态 修复 和 土地 整治 策略 
提供 科学 借鉴 。 但 是 本 研究 是 在 30 m 分 辩 率 遥感 
数据 的 基础 上 完成 的 地 类 解 译 工作 ,研究 过 程 存 在 
一 定 的 主观 判断 ,不 利于 对 土地 利用 变化 和 生境 质 
量 开 展 精 准 评估 ,因此 ,今后 的 研究 应 加 强 多 尺度 
综合 性 研究 ,提升 解 译 精度 ,为 强化 结果 精度 提供 
更 加 坚实 的 数据 基础 。 


4 结论 
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时 间 序 列 下 的 土地 利用 信息 ,结合 InVEST 模 型 分 析 
1990—2020 年 拉萨 河流 域 土地 利用 和 生境 质量 时 
空 变化 特征 ,得 到 如 下 结论 : 

(1) 针对 高 寒 高 海拔 地 区 的 土地 利用 分 类 问 
题 ,进行 遥感 数据 可 用 性 分 析 是 有 必要 的 ,结果 表 
明 :筛选 植被 生长 季 , 用 前 后 年 份 的 影像 进行 补充 
和 替换 缺失 的 数据 ,并 去 云 后 生成 的 年 度 合 成 影 
像 ,提高 了 地 类 分 类 效率 和 精度 ;此 外 ,建立 了 一 种 
结合 Landsat 影 像 光谱 特征 纹理 特征 和 地 形 特 征 的 
最 优 RF 分 类 模型 ,实现 了 对 拉萨 河流 域 土地 利用 
信息 的 分 类 ,OA、Kappa 系数 和 Fl-score 分 别 在 
88.16% .0.84 和 0.70 以 上 。 

(2) 拉萨 河流 域 草地 分 布 最 广 ,多 年 平均 占 比 
2) 56.74% ,其 次 是 未 利用 地 , 约 占 28.33%。20 世纪 
90 年 代 以 来 ,拉萨 河流 域 土地 利用 变化 特征 整体 表 
现 为 "五 增 三 减 ” ,草地 建设 用 地 、 林 地 湿地 和 水 
域 面积 增加 ,耕地 、 冰 川 及 永久 积 雪 和 未 利用 地 面 

只 减少 ,其 中 ,建设 用 地 的 增长 幅度 最 大 ,达到 
288.35% , 而 耕地 缩减 最 多 , 达 50.18% , 且 耕 地 减少 
主要 是 受到 城镇 范围 的 挤占 ;未 利用 地 面积 减少 最 
为 显著 ,其 转 出 部 分 主要 转化 为 草地 ,表明 气候 变 
化 和 生态 修复 工程 对 拉萨 河流 域 生 态 恢 复 均 起 到 
了 推动 作用 。 

(3) 1990 一 2020 年 拉萨 河流 域 生境 质量 从 空间 
分 布 上 来 看 ,高 质量 生境 主要 分 布 流域 下 游 .西南 
侧 以 及 源头 区 域 ,以 草地 ,水体 .冰川 及 永久 积 雪 和 
湿地 等 土地 利用 类 型 为 主 ; 而 拉萨 市 市 辖区 A HC 
区 、 林 周 县 . 曲 水 县 和 墨 竹 工 卡 县 县 城 周 边 ,以 及 流 
域 中 上 游 的 未 利用 地 等 区 域 生境 质量 较 差 。 从 数 
量变 化 上 来 看 ,研究 期 内 流域 整体 生境 质量 呈 " 先 
减 后 增 再 减 然 后 趋 缓 "变化 ,生境 质量 指数 从 0.53 
提升 到 0.57。 

(4) 研究 识别 了 1990 一 2020 年 土地 利用 变化 
和 生境 质量 变化 的 主要 原因 ,人 口 增加 和 快速 城市 
化 加 快 了 城镇 范围 扩张 和 耕地 的 减少 ,但 气候 变化 
和 生态 工程 对 拉萨 河流 域 的 生境 质量 改善 均 起 到 
了 积极 作用 。 
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Spatial-temporal evolution characteristics of land use change and habitat 
quality in the Lhasa River Basin over the past three decades 
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(1. College of Ecology and Environment, Southwest Forestry University, Kunming 650224, Yunnan, China; 
2. College of Forestry, Southwest Forestry University, Kunming 650233, Yunnan, China) 


Abstract: This study used the Lhasa River Basin, which is a typical plateau basin to clarify the spatiotemporal 
characteristics of land use change and its impact on a plateau basin habitat quality. The google earth engine 
(GEE) platform and random forest classification algorithm were used to interpret the land use information of the 
Lhasa River Basin from 1990 to 2020, and to analyze the spatiotemporal dynamic changes of land use in the past 
30 years. The INVEST model was integrated to evaluate the habitat quality of the basin, and to explore the 
characteristics of its habitat quality changes. The results showed that the efficiency and accuracy of land use 
classification improved by selecting the vegetation growing season, supplementing and replacing missing data 
with images obtained in the previous and following years, and generating annual composite cloud free images for 
land use classification in the alpine and high-altitude areas. An optimal RF classification model combining the 
spectral, textural, and topographical features of Landsat images was established to classify land use information 
in the plateau basin. The observed OA, Kappa coefficient, and the Fl- score were 88.16%, 0.84, and 0.70, 
respectively. Since the 1990, the features of the land use change in the Lhasa River Basin have been characterized 
by “five increases and three decreases,” with an increase being observed in the areas under grassland, 
construction, forest land, wetland, and water, while the areas under cultivated land, glaciers, and permanent snow 
or unused land having a decrease. Among them, the growth rate of construction land had the largest increase of 
288.35% , while the largest decrease was detected in cultivated land with 50.18%, which could mainly be 
attributed to overpopulation in the urban areas. Most significant decrease occurred in the unused land area, and 
the changed area was mainly converted with grassland. During the study period, high quality habitats were 
mainly distributed in the downstream areas, southwest, and the main source areas of the basin, included the 
grassland, water, glaciers, permanent snow, and wetland. Poor qualities of habitats were observed in the Lhasa 
municipal district, Dazi District, Linzhou County, Qushui County, and Mozhugongka County, as well as in the 
unutilized land in the middle and upper reaches of the basin. The overall habitat quality of the watershed showed 
an alternating transition pattern of initial decrease, then increase, followed by a decrease, and a slowed down 
transition, while the habitat quality index increased from 0.53 to 0.57. Population growth and rapid urbanization 
accelerated the expansion of urban areas and the reduction of arable land, while climate change and ecological 
engineering played a positive role in improving the quality of habitat in Lhasa River Basin. This study provides a 
scientific basis and reference for formulating ecological restoration strategies, sustainable land use, and habitat 
quality improvement in a plateau area. 
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